Dragging D mesons by hot hadrons 
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We evaluate the drag and diffusion coefficients of a hot hadronic medium consisting of pions, 
nucleons, kaons and eta using open charm mesons as a probe. The interaction of the probe with 
the hadronic matter has been treated in the framework of effective field theory. It is observed that 
the magnitude of both the transport coefficients are significant, indicating substantial amount of 
interaction of the heavy mesons with the thermal hadronic system. The results may have noticeable 
impact on the experimental observables like the suppression of single electron spectra originating 
form the decays of heavy mesons in nuclear collisions at relativistic energies. 
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Nuclear collisions at Relativistic Heavy Ion Col- 
lider (RHIC) and Large Hadron Collider(LHC) en- 
ergies are aimed at creating a phase where the 
bulk properties of the matter are governed by a de- 
confined state of (light) quarks and gluons known 
as Quark Gluon Plasma (QGP). The study of the 
transport properties of QGP is a field of great con- 
temporary interest and the heavy flavors, namely, 
charm and bottom quarks, play a crucial role in 
such studies. The weakly interacting picture of the 
QGP stems from the perception of asymptotic free- 
dom of QCD at high temperatures and densities. 
However, the experimental data from RHIC, espe- 
cially the measured elliptic flow indicate that the 
matter produced in Au+Au collisions exhibit prop- 
erties which are more like a strongly interacting liq- 
uid than a weakly interacting gas. The magnitude of 
the transport coefficients can be used to understand 
the strength of the interaction within the QGP. For 
example, the shear viscosity or the internal friction 
of the fluid symbolizes the ability to transfer mo- 
mentum over a distance of about a mean free path. 
Therefore, in a system where the constituents inter- 
act strongly the transfer of momentum is performed 
easily - resulting in lower values of r/. Consequently 
such a system may be characterized by a small value 
of rj/s where s is the entropy density. On the other 
hand, for a weakly interacting system the momen- 
tum transfer between the constituents become stren- 
uous which gives rise to a large rj. The importance 
of viscosity also lies in the fact that it damps out 
the variation in the velocity and makes the fluid flow 
laminar. A very small viscosity (large Reynold num- 
ber) may make the flow turbulent. A lower bound 
on the value of rj/s has recently been found using 
AdS/CFT Q. 

The interaction of heavy quarks with the QGP 
can be used to estimate the value of the transport 
coefficients. This has recently been performed by ex- 
perimentally measuring the nuclear suppression fac- 
tor (Raa) 0| and the elliptic flow (V2) [3] for the 
single electron spectra originating from the semi- 



leptonic decays of the heavy mesons which are pro- 
duced from heavy quark fragmentation. Several the- 
oretical attempts [4H15| have been made to explain 
Raa and vi , where the role of hadronic matter has 
been ignored. However, to make the characteriza- 
tion of QGP reliable the role of the hadronic phase 
should be taken into consideration and its contribu- 
tion must be subtracted out from the observables. 
Although a large amount of work has been done on 
the diffusion of heavy quarks in QGP, the diffusion of 
heavy mesons in hadronic matter has received much 
less attention so far. Recently the diffusion coeffi- 
cient of D meson has been evaluated using heavy 
meson chiral perturbation theory [l6| and also by 
using the empirical elastic scattering amplitudes [17| 
of D mesons with thermal hadrons. 

The present work addresses the relevance of the 
hadronic sector to some of these issues. The drag 
and diffusion coefficients for the hadronic phase have 
been evaluated and their importance to the experi- 
mental observables in heavy ion collisions have been 
discussed. In particular, we consider the interaction 
of a D meson with a thermal hadronic system com- 
posed of pions, nucleons, kaons and 77 in a tempera- 
ture domain relevant for heavy ion phenomenology. 
It is expected that the relaxation time for heavy 
mesons are larger than the corresponding quanti- 
ties for light hadrons. The abundance of the heavy 
mesons will be low for the temperature (T) range 
under consideration (T=100-180 MeV), as a result 
they do not decide the bulk properties of the matter. 
The thermal production of charm mesons can be ig- 
nored for the range of temperature mentioned above. 
Therefore, the drag (7) and diffusion (Bo) coeffi- 
cients of the heavy mesons can be evaluated by using 
its elastic interaction with the thermal hadrons. For 
the process, D(p) + h(q) — > D(p r ) + h(q') (h stands 
for pion, nucleon, kaon and eta), the drag 7 can be 
calculated by using the following expression (l8j : 
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go being the statistical degeneracy of the D meson 
propagating in the medium. The above expression 
indicates that the drag coefficient is the measure of 
the thermal average of the square of the invariant 
amplitude | M | 2 weighted by the momentum trans- 
fer, p—p'. The factor f(q) denotes the thermal phase 
space for the particle in the medium. 

Similarly the diffusion coefficient Bq can be de- 
fined as: 
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With an appropriate choice of T(pf) both the drag 
and diffusion co-efficients can be evaluated from the 
following expression: 
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where A(x, y,z) = i 
is the triangular function and | M | 2 in the present 
case corresponds to the scattering of O-mesons from 
the light mesons n, K, r\ and nucleons. In a hot 
pion gas this was obtained by Fuchs et al. 
using experimental information on Z?-meson reso- 
nances. Here we evaluate these amplitudes using a 
covariant formulation of chiral perturbation theory 
in which the leading term allows for D scattering 
via D* meson exchanges in addition to a contact 
interaction. In nuclear matter the DN scattering 
amplitudes have been obtained in a coupled chan- 
nel Bethe-Salpeter approach where A c and S c ap- 
pear as dynamically generated [20] . In this work 
we have obtained the DN scattering amplitudes 
proceeding via A c and S c exchanges using the La- 
grangian of Ref. [2l[. The generic Feynman dia- 
grams for the elastic processes are depicted in Fig. [TJ 
We have included form factors in each of the inter- 
action vertices to take into account the finite size of 
the hadrons. For the t and s-channel diagrams the 
form factors are taken as [2l[ F t — A 2 / (A 2 + q 2 ) and 
F s = A 2 / (A 2 +pl ) respectively, where q is the three 
momentum transfer and pi is the initial three mo- 
mentum of the light mesons (pion, kaon and eta) or 



nucleon. In the four point (contact) vertices a form 

factor, F 4 = (A 2 /(A 2 + q 2 )) 2 with q = 5p 2 /3 has 
been introduced [22J. We have taken A = 1 GeV. 
The interaction Lagrangian as well as the | M | 2 for 
various processes are detailed in the appendix. 





FIG. 1: Feynman diagrams for the scattering of D 
mesons with hadrons in the medium. Here M stands 
for mesons (pion, kaon and eta) and P* and B denote 
charmed vector mesons and baryonic resonances respec- 
tively. 



In Fig. [^variation of the drag coefficient with tem- 
perature has been depicted for D-mesons. We have 
observed that the D — ir meson interaction plays the 
most dominant role in the drag coefficient primarily 
because of the larger phase space density of the pi- 
ons. However, at higher T the contribution from the 
nucleons become significant. As mentioned before, 
7 is the thermal average of the square of the invari- 
ant amplitude weighted by the momentum transfer. 
Therefore, as the temperature of the thermal bath 
increases the hadrons move faster and gain the abil- 
ity to transfer larger momentum during their inter- 
action with the D mesons - resulting in the increase 
of the drag coefficient. This trend is clearly observed 
in Fig. [5J It may be mentioned here that the drag 
increases with T when the system behaves like a gas. 
In case of a liquid the drag may decrease with tem- 
perature (except for very few cases) since a substan- 
tial part of the thermal energy goes into making the 
attraction between the interacting particles weaker. 
This allows them to move more freely resulting in a 
smaller drag force. Therefore, the variation of the 
drag with T may be used to characterize the nature 
of interaction of the fluid. The large value of the 
drag coefficient indicates that the interaction of the 
D meson with the thermal medium is quite signifi- 
cant so that the the D meson may get thermalized 
in the system and flows with the bulk matter. This 
may be examined by analyzing the transverse mo- 
mentum spectra of the D mesons produced in heavy 
ion collisions [23 |, We find that the typical value 
of the relaxation time (~ 7" 1 ) is about 4-5 fm/c. 
Therefore, if the life time of the hadronic phase is 
more than this time scale then the D meson may 
get thermalized indeed. 

In the same way it may be argued that the diffu- 
sion coefficient involves the square of the momentum 
transfer - which should also increase with T as seen 
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FIG. 2: The variation of drag coefficients with tempera- 
ture due to the interaction of the D with thermal pions, 
nucleons, kaons and eta. 
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FIG. 3: Variation of diffusion co-efficient as a function 
of temperature. The solid line indicates the variation of 
the diffusion coefficient with temperature obtained from 
Eqs. [3] and [4] The dashed line stands for the diffusion 
coefficient obtained from the Einstein relation (Eq. [5}. 



in Fig. |31 The dominant contribution comes from 
the interaction of D mesons with pions. The drag 
and the diffusion coefficients are related through the 
Einstein relation as: 

B = M D1 T. (5) 

where Md is the mass of the D-meson. The tem- 
perature variation of the diffusion coefficient ob- 
tained from Eq. [5] is depicted in Fig. [3] (dashed 
line). The difference between the results obtained 
from Eq. U and Einstein's relation is about 10% at 
T = 180 MeV. This small difference illustrates the 
validity of the Einstein relation in the current situ- 
ation. The value of the spatial diffusion coefficient, 
D x may be expressed in terms of drag coefficient as 
D x = T/{M Dl ). The value of D x at T = 180 MeV 
is ~ 2/(27rT) i.e. 2 times larger than the thermal 



wave length, A = 1/(2%T) which is well within the 
quantum bound. 

The magnitude of the energy dissipation of the D 
meson in the system may be estimated by using the 
relation 



The magnitude of 7 obtained in the present calcula- 
tion indicates a substantial loss of energy of the D 
meson in the medium, which might have observable 
effects on quantities such as the nuclear suppression 
factor of single electrons originating from the decays 
of heavy mesons. 

Here it is necessary to point out that though chiral 
perturbation theory provides a consistent framework 
for performing perturbative calculations of strong in- 
teraction processes such as D — h scattering in this 
case, it is limited by the abundance of coupling con- 
stants appearing in the Lagrangian which have to 
be determined from experimental data (24|. In the 
case at hand, the experimental error in the D* decay 
width leads to an uncertainty (11%) in the D*Dn 
coupling, g [25| . which results in a significant vari- 
ation in the value of the drag diffusion coefficients. 
Interestingly, the lower bound in this value leads to a 
drag coefficient which agrees reasonably well at high 
temperatures with that obtained by He et al [TtJ us- 
ing empirical elastic scattering amplitudes. In addi- 
tion to this, chiral symmetry breaking effects in the 
pseudoscalar decay constants could also contribute 
to the uncertainty. However, since in our case the 
pion contribution dominates maximally, this effect 
will be insignificant. 

To summarize, in this work we have evaluated the 
drag and diffusion coefficients of open charm mesons 
propagating in a hadronic background composed of 
pions, kaons, nucleons and eta. We observe that the 
values of both the transport coefficients increases 
with temperature and the dominant contributions 
come from the pions in the medium. However, at 
higher T the contributions from heavier hadrons be- 
come significant. The magnitude of the drag coef- 
ficient of the D meson in the hadronic medium re- 
veals that while evaluating the nuclear suppression 
for the single leptons originating from the decays of 
D mesons the hadronic contributions should be in- 
cluded. Lattice QCD calculations [26[ indicate that 
at low baryonic chemical potential and high temper- 
ature domain there is no phase transition between 
hadronic matter and QGP - it is a cross over, which 
means that the hadronic matter can make a con- 
tinuous transition to QGP in this region of phase 
diagram. Therefore, the transport coefficients eval- 
uated for the hadronic matter with zero baryonic 
chemical potential may have vital effects from the 
quark gluon plasma. 
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I. APPENDIX 



In this appendix we provide the interaction La- 
grangian and matrix elements for scattering of D 
mesons from the light mesons (tt, K, g) and nucleons 
discussed in this work. 

The leading order chiral Lagrangian describing the 
interaction of Goldstone bosons with the heavy-light 
pseudoscalar (P) and vector (P* ) mesons is given by 



where P 



+ig(P*u^P Jf - Pu^P l 

{D°,D+,D+) and 



(P*»P?) 
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(-D* , D* + , D*+) and (...) denotes trace in flavour 



space 

V»P a = d^P a 



The covariant derivatives are defined as 



p 6 it° and v pt = a^pt + r^p. 
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with a, b the Sf7 (3) flavour indices. The value of the 
heavy-light pseudoscalar-vector coupling constant 
g = 1177 ± 137 MeV is obtained by reproducing 
the experimental D* —> Dir decay width of ~ 
65±15 keV with the above interaction. The vector 
and axial-vector currents are respectively given by 
= ^(u^d^u + udfjrf) and — i(ytdpU — ud^v)) 
where u = exp(^-). The unitary matrix $ 
collects the Goldstone boson fields and is given by 



V2 



V 



7T 

K+ 



it" I V 

K° 



K° 



To lowest 



2n 

V6 



order in $ the vector and axial-vector currents are 
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The above interaction allows elastic scattering of 
the D meson with the tt, K and g fields through a 
heavy-light vector meson exchange in addition to a 
contact interaction as shown in Fig. ([1]). The form 
of the contact interaction obtained in the covari- 
ant formulation of chiral perturbation theory used 
here coincides with that in |27| . The invariant am- 
plitudes for D + clastic scattering from hadrons, h 
(D+(pi) + h(p 2 ) -> D+(p 3 ) + h(p 4 )) have been ob- 
tained as follows: 
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In the numerical calculations we have used the 
physical masses for the particles involved. For the 
heavy-light mesons we have taken mo = 1867 MeV, 
m D ' = 2008 MeV and m D , = 1969 MeV. 

We next discuss scattering of the D meson 
with nucleons which proceeds via exchange of the 
charmed baryons E c and A c according to the La- 
grangian plj . 
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where D = (D°,D+) and D = (D°,D-) T . The am- 
plitudes for D + mesons elastically scattering from 
nucleons (and anti-nucleons) via A c exchange are ob- 
tained as, 
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The coupling constants obtained using SU(A) sym- 
metry are given by [U| = 7.18 GeV" 1 and 
=2.01 GeV- 1 . 
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and 



The various invariant amplitudes can be expressed 
in terms of the Mandelstam variables using the re- 

-i ,. s — m 2 n —mi m 2 n +m h — t 

lations: pi ■ p 2 = § -, Px • P3 = 2 

2 , 2 

n+ m h 



and pi ■ Pi 



where is the mass of D 



meson and mu that of the light hadrons. 
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